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We study multipartite correlations and non-locality in an isotropic Ising ring under transverse 
magnetic field at both zero and finite temperature. We highlight parity-induced differences between 
the multipartite Bell-like functions used in order to quantify the degree of non-locality within a ring 
state and reveal a mechanism for the passive protection of multipartite quantum correlations against 
thermal spoiling effects that is clearly related to the macroscopic properties of the ring model. 
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I. INTRODUCTION 

Correlated systems are fundamental in developing and 
exploiting the possibilities arising from quantum mechan- 
ics. A premier role in such an endeavor is played by 
chains of interacting quantum spins [l|, in light of their 
wide applicability to many physical settings and their 
ability to simulate the behavior of strongly correlated 
systems in condensed matter 0, 0] . They are a key in- 
gredient in many distributed quantum computation and 
communication protocols If we are truly to benefit 
from the potential of such resources, a thorough under- 
standing of the nature of such correlations is essential. In 
this respect, it is particularly important to gather infor- 
mation on the distribution and amount of genuine mul- 
tipartite quantum correlations in one of such systems so 
as to be able to tailor specific protocols to the available 
structure of entanglement among the parties of a chain 
or lattice. 

Entanglement and non-locality are the key figures of 
merit for non-classical correlations. Already a good un- 
derstanding exists for simple bipartite entangled systems 
and there is much research in the direction of multipar- 
tite entangled states [|. Equally, non-locality is truly at 
the forefront of current research as the premier manifes- 
tation of true deviation of the behavior of a system from 
classicality Non- locality embodies a stronger type 

of correlation than entanglement and the development of 
tools to detect, quantify and characterize it is moving 
swiftly in the field. Of course one should be mindful of 
how any such system can be experimentally generated 
and utilized and of the fact that unavoidable interac- 
tions with a surrounding environment result in the loss 
of such strong quantum correlations. Methods to con- 
trol the affect of these unwanted interactions (resulting 
in either losses or pure decoherence) have been proposed, 
ranging from bang-bang tech niq ues @ to the active use 
of the quantum Zeno effect flCfl. All such methods, to 
the best of our knowledge, require either heavy infiu- 
ences on the system at hand or expensive additional re- 
sources for entanglement protection. The generation of 
non-classical multipartite correlated states is already a 
formidable task: if loss/decoherence-tolerant methods for 
their production have to be in order, this would neces- 



sarily imply a frequently prohibitively escalating exper- 
imental cost. It should thus be evident that a scenario 
where highly quantum correlated states are produced as 
the ground state of a given many-body interacting model 
and naturally maintained while in the proximity of envi- 
ronmental effects, is endowed with extreme physical ap- 
peal. 

In this paper we present a study of the Ising model 
of spin-1/2 particles (qubits) in a transverse magnetic 
field. This model has attracted wide attention due to its 
exactly solvable nature and the attractive natural prop- 
erties of its ground states [lTI - [l3| . We go significantly 
further than the available literature on this model so 
far by addressing the case of genuine multipartite quan- 
tum correlations and their nature. We examine the zero 
temperature dynamics of entanglement and non-locality 
highlighting by means of an efficient hybrid approach 
for the determination of multi-site correlation functions. 
The finite temperature regime is then addressed and we 
show the ability to control and preserve the non-locality 
present in the system with minimal additional effort. We 
show the key ingredient in the model is the presence and 
manipulation of the magnetic field. 

The remainder of the paper is organized as follows. 
In Sec. |lT]we introduce the physical model for spin-spin 
coupling considered in this work and sketch the technique 
required for its full diagonalization. Sec. lIIII discusses the 
technical tools used in our investigation for the quantifi- 
cation of multipartite entanglement and quantum non- 
locality. Sec. I IV I presents the dynamics of multipartite 
correlations at zero temperature while Sec. |V] we con- 
sider how entanglement and non-locality are affected in 
the finite temperature case. Finally Sec. IVII summarizes 
our findings. 



II. THE MODEL 

We begin our study by discussing features and proper- 
ties of an isotropic Ising model in a transverse magnetic 
field. For N coupled spins the Hamiltonian model reads 

N N 

u^-jY.K^K+i + bY.K (1) 

rt=l n—1 
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with J the inter-spin couphng strength and B the global 
magnetic field for each qubit and a\ (i = x,y,z) the 
i Pauli spin operators for spin n — 1,..,N. In order 
to strip our analysis from unnecessary complications, we 
work with a dimensionless version of the Hamiltonian 
defined by assuming J = J8 and B = B£ with £ a 
common order-of- magnitude factor and J, B are dimen- 
sionless parameters. We thus consider the dimensionless 
Hamiltonian 

H^^^-jY.^n®K+i+BY,K (2) 

n— 1 n— 1 

We can readily determine the thermal state associated 
with this Hamiltonian using 

e{T)^e-^^/Z (3) 

where Z is the partition function and /3 — h/ks- Clearly 
each spin only interacts with its nearest neighbors and 
we shall assume a cyclic boundary condition such that 
the last spin interacts with the first, thus forming a spin- 
ring. It is also worth noting that we have taken the mag- 
netic field to act orthogonal to the interaction. This par- 
ticular model is very attractive due to its exactly solv- 
able nature. The solution is presented in [ll|, [l^ [T6| . 
and so for completeness we simply comment on the 
methods employed. In order to diagonalize the Hamil- 
tonian we use the Jordan- Wigner transformation jl7| . 
Firstly, we transform the Hamiltonian using a new vari- 
able cr^ = 1/2 (cr^ ± ial^). We can then easily move into 
the fermionic picture by introducing a further transfor- 
mation 

4 = (cn)^ = ®"jZl {c„, 4} = 5„„, (4) 

with 5nm the Kronecker delta, which is 1 for n = m and 
otherwise. These fermionic variables then allow the 
Hamiltonian to be written in the form 

N N ^ 

iJ=-J^(ct-c„)(cl+i-|-c„+i)+2i?^(ctc„--). (5) 

n— 1 n— 1 

Using a Fourier transform we pass to the momentum rep- 
resentation. The resulting Hamiltonian can then be di- 
agonalized by means of a Bogoliubov transformation in- 
troducing new fermionic operators {6^,6]^}. The latter 
are related to the Fourier-transformed fermionic opera- 
tors c±k through bk — cos(i?i./2)cfc — j sin(i?fc/2)cLj, and 
— cos(i?fc/2)c[.-|-« sin(i?j./2)c_fc. This brings it into the 
free-fermion form 

Hif = J2^kblbk-J2^k, (6) 

fe A; 

where = J'^ + B"^ — 2JB cos 0fe , tani?fe = {—B + 
J cos (j)k)/ {J sm(t)k) and (pk — TT{2k + 1)/N with k — 
—N/2,..,N/2 — 1 in the subsector with an even num- 
ber of fermions and a slightly different definition for the 



odd-number case. However, in what follows we will be 
mostly concerned with the ground state of the Hamilto- 
nian, which lives in the even-number sector (the appro- 
priate expression will be used when leaving the ground 
state). The constant term in Eq. (|6]) is related toj^k ^k- 
The ground state of the system is then the state with no 
population in any of the decoupled bk modes, i.e. the 
state solves the set of conditions bk IgSjy) = simultane- 
ously, for any allowed value of fc. A standard calculation 
(see for details Ref. [l^) then leads to the ground state 

|gsw) =(8)[cos^|00),__, + zsin^|ll),__,] (7) 
fc 

with |0)^^ {\^) ^f.) the state with (1) fermions with 
momentum (pk- This approach is particularly useful 
when a ring with a large number of particles is consid- 
ered. In this case, the direct diagonalization of Eq. 
in fact, becomes computationally demanding. Regard- 
less of the dimension of the ring at hand, the mag- 
netic field proves to be key to the behavior of entangle- 
ment and non- locality within our system. In fact, when 
no magnetic field is present, from Eq. (|6]) we see that 
Ck = J V/c, so that the ground state has energy —NJ, 
which is the same for (^^=1 !+)« ($n=i \~)n '^'^^^ 
|±)„ = (|0)„ ± |l)„)/\/2 being the eigenstates of a^. In 
these expressions, {|0)„ , |1)„} is the computational logic 
basis of spin n — 1, .., N . As Stelmachovic and Buzek 
proved in Ref. the proper ground state of the ring 
in these conditions can be taken as a Greenberger-Horne- 
Zeilinger state [14] in the basis built from the tensor prod- 
uct of eigenstates of a^. That is 

|GHZ^) = (|+ + ..+) + I- -..-))/V2, (8) 

which, when written in the ®n=i^n eigenbasis, always 
has an even number of spins in their |1) state. In what 
follows, the energy of the ground state on a system of N 
spins is indicated as 

Ajv = -^efc. (9) 

k 



III. TOOLS FOR STUDYING MULTIPARTITE 
ENTANGLEMENT AND NON-LOCALITY 

In order to examine the dynamics of non-locality and 
entanglement, we need appropriate instruments to quan- 
tify them. Bipartite entanglement can be quantified by 
a number of well-defined measures. However for mixed, 
multipartite states no such measure exists [3]. Here we 
we shall employ a suitable generalization of negativity, 
an entanglement measure based on the well-known Peres- 
Horodecki criterion [l^ , which is necessary and sufficient 
for 2 X 2, 2 X 3 and oo x c» systems . , together with an 
extension of Wootters concurrence [20| to the case of a 
multipartite register of an even number of elements. 
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For a two-qubit system, negativity is defined as 

AA2 = -2max[0,Anog], (10) 

where Ancg is the single negative eigenvalue of the par- 
tially transposed version of the density matrix describing 
the state of the qubits. This measure can be extended 
to tripartite systems by simply considering the partial 
transposition of the density matrices associated with ap- 
propriate bipartitions of a given system (2]| . Let us con- 
sider three spins. A, B and C. The entanglement of spin 
A with spins B and C is 

^A(BC) = -2max[0, ^ A^^lff ]. (11) 
j 

with the negative eigenvalues of the partially 

transposed density matrix describing the bipartition A — 
{B,C). We can then determine the genuine tripartite 
entanglement content of the tripartite state by taking a 
geometric average of all the bipartitions, thus ensuring 
that the corresponding measure remains a true entangle- 
ment monotone. Therefore 

^f3 = [^fA{BC)^fB{AC)^fciAB)]^ ■ (12) 

Unfortunately, the extension of this measure to larger 
systems is not as straightforward and one has to resort 
to other tools for the quantitative analysis of multipar- 
tite entanglement. A rather handy one is embodied by 
iV-concurrence ^2§\ , which is the direct generalization of 
bipartite concurrence to the case of an even number of 
qubits. That is, one has to consider the eigenvalues rjj of 
the A'^-spin flipped matrix 

P«=i^Pp*«=i'^P (13) 

with p the density matrix describing the state of the sys- 
tem and p* its complex conjugate. The A^-concurrence 

NC is then given by NC — max[0, v^i~X]^2 with 
771 > rjj \f j. This measure is appealing due to its compu- 
tationally handy form and its sensitivity to global entan- 
glement. In fact, it is exactly zero if any qubit is separa- 
ble from the rest of a system. Its non-nullity, though, is 
only a sufficient condition for multipartite entanglement 
as there are multipartite entangled states with vanishing 
A^-concurrence, such as the N-qubit H^-states. According 
to the definition given by Wong and Christensen j223 , in 
pure states N-concurrence quantifies the overlap between 
a given state and the one obtained upon application of 
an anti-unitary time-reversal operation. For qubit states, 
such operation is embodied by a collective spin-flip oper- 
ation. A-concurrence is thus a natural choice when the 
state to study is the eigenstate of a Hamiltonian respect- 
ing time-reversal symmetry, such as those that can be 
written as the sum of tensor products of an even number 
of non-identity Pauli matrices [2^. The Ising Hamilto- 
nian satisfies such condition when no magnetic field is 
applied, thus justifying our choice and our efforts here 



can well be seen as the quantification of the robustness 
of the ground state of a Ising-like model to the perturba- 
tion induced by the magnetic field. 

There is no formal way to quantify the degree of 
non-locality of a given state. Nevertheless, the degree 
of violation of Bell-type inequalities can be taken as a 
an intuitive means of ascertaining the amount of non- 
locality (2^ . A more pragmatic viewpoint to quantify 
non-locality, on the other hand, is to consider the amount 
of noise required by a given state to stop violating a Bell- 
like inequality. Moreover, in the multipartite setting, a 
few ambiguities still remain in deciding which type of in- 
equality to utilize. In fact, as very clearly pointed out 
by Cereceda |25| by reprising a point originally made by 
Svetlichny '26|, n-particle entanglement is not physically 
equivalent to n-particle non-locality. In this respect, the 
violation of an n-particle Bell-like inequality of some sort 
by an n-particle entangled state is not enough, per se, 
to prove genuine multipartite non-locality. When one 
moves into the multipartite setting Bell-type inequalities 
exist that probe various different types of correlation. 
Technically, this issue has produced considerable theo- 
retical effort directed towards the formulation of appro- 
priate Bell-like inequalities able to capture the genuine 
multipartite non-local nature of a given state. Notice- 
ably, Svetlichny has derived an inequality for the tripar- 
tite case that, while obeyed by models assuming two- 
particle non-locality, is violated by quantum mechanical 
states that are genuinely three-particle entangled [26] . As 
discussed by Cereceda [1^, a Svetlichny inequality (SI) 
is a righteous Bell inequality for the tripartite case and 
emerges as a valuable tool for the unambiguous asser- 
tion of the existence of genuine tripartite entanglement 
and tripartite non-locality in any three-particle state. To 
this task, the use of a standard Mermin's inequality [271 is 
not sufficient: on the contrary, quantum correlations that 
violate SI are strong enough to maximally violate Mer- 
min's inequality as well. The inequality by Svetlichny 
has been independently extended to the n-partite sce- 
nario in Refs. [Mllil- In particular. Coffins et al. [2!] 
have provided an iterative way to construct multipartite 
Bell inequalities for dichotomic observables. We shall call 
Oj and Oj the two outcomes of a local measurement per- 
formed over one of the particles. It should be noted, the 
explicit form of the observable can be dependent on the 
state to be tested. For instance, measuring along the in- 
correct axis for a particular state can lead to a failure to 
see any violation. In the proceeding section we will fix the 
appropriate observable for our model. Then, by setting 
nil = oi (Ml = Oi), one constructs the polynomials 

TO„ = ■^m„_i(o„ -I- On) + \Mn-l{On - On), 
M„ = -M„_i(o„ -I- On) + -jmn-l{On - 0„). 

Quantum mechanically, we shall interpret the polynomi- 
als rn„ and Af„ as sums of n-particle correlation functions 
for measurements having outcomes Oj and Oj. We can 
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FIG. 1: (Color online). Behavior of entanglement against B 
and J for (a) bipartite spin ring and (b) tripartite spin ring. 
For small values of B the ground state is highly entangled. 
As the field strength is increased we see the entanglement 
asymptotically approach zero. 



then define the generalized SvetUchny polynomials 



Sn — 



nin {n even) 

(m„ + Af„)/2 (n odd) 



(15) 



The bound imposed to 5„ by local hidden- variable mod- 
els is 1, while quantum mechanically an n- qubit GHZ 
state achieves the maximum value of V2"~^ and v2^^ 
for an even and odd number of particles, respectively. 
Armed with such tools, we now tackle the main point of 
our study. 



IV. BEHAVIOR OF NON CLASSICAL 
CORRELATIONS: ZERO-TEMPERATURE CASE 

We begin our analysis considering the dynamics of the 
correlations in the simplest case of just two coupled spins. 
By assuming that both the spins are at zero temperature, 
the ground state of Eq. 1^ has the form 



?S2) 



-(i3 + A2/2)|00)i2 + J|ll)i2 



J2 + (B - 



(16) 



with A2 = -2\/B2 + J2 > B. Clearly, for B ^ 0, 
IgSj) (|00) + |ll))/\/2 = IGHZ2), in agreement with 
the expectations for the ground state of the model at 
hand. The calculation of negativity is therefore straight- 
forward and leads to 



^^2 



J 



+ J2 ^ 



(B > 0). 



(17) 



The behavior against J and B is shown in Fig. [T] (a). 
Clearly, as S — >■ the entanglement achieves a full cbit, 
while as we increase the amplitude of the magnetic field, 
the amount of zero-temperature entanglement decreases, 
vanishing asymptotically. This behavior is in agreement 
with what has been shown in Refs. [lll.[T2j. where concur- 
rence was considered. We can extend our study to the 
case of tripartite entanglement. The zero-temperature 



tripartite negativity of the ground state of a three-spin 
ring can be calculated using Eq. ([T^ to be 



^f3 



1 



3J2 



4S-2J 



■iV B-2 - BJ + J-2 -BJ+J^ 



+ 4. (18) 



Fig. [T] (b) shows its behavior. Analogously to the en- 
tanglement of the two-spin ring, we see that as the mag- 
netic field is increased there is an asymptotic decay to 
zero. While a.t B the state is the maximally en- 
tangled GHZ state, interestingly as we increase the mag- 
netic field strength we find the reduced state of any two 
spins becomes entangled, thus stating the departure of 
the ground state of the system from the GHZ class of 
entangled states. As discussed in Sec. Illli the case of 
an arbitrarily sized ring is not easily managed because 
of the lack of available measures for genuine multipar- 
tite entanglement. Yet, partial information on the dis- 
tribution of quantum correlations within an even-parity 
length ring can be effectively gathered by means of the 
TV-concurrence. We have thus evaluated the NC for the 
significant cases of iV = 4 and 6 and compared it with the 
entanglement measured by Wootters concurrence evalu- 
ated for a two-spin chain. The results are shown in Fig. [5] 
against J and B and reveal an increasing fragility of mul- 
tipartite entanglement as N grows. As it should be ev- 
ident, at non-zero J entanglement only asymptotically 
decreases as B grows showing that, although the ground 
state of the ring progressively deviates from a GHZ state, 
a degree of multipartite entanglement is nevertheless kept 
by the spins. 

However, the inability to discuss the non-classical fea- 
tures of correlations in the odd-parity scenario and the 
lack of a "homogeneous" way to quantify entanglement 
(two different measures had to be used in the above) is 
certainly an unsatisfactory feature. Therefore, in order 
to provide a much more self-consistent characterization of 
the correlation properties within the ground states of the 
Ising ring, we move to the study of non-locality. In such a 
two-spin setting, the well-known Closer-Horne-Shimony- 
Holt (CHSH) inequality [23| can be used, which is en- 
compassed in the formulation of the general multipartite 
problem given in Sec. IIIII In order to run the CHSH in- 
equality test, we thus project the state of the two-spin 
system onto the cigenstates of the local observable 



cos f 
i sin( 



—I smf 
— cos 6 



(z-1,2) 



(19) 



and calculate the correlation function among the various 
combinations of outcomes. This is equivalent to 



C{9i,e2) = {gs^\Ri{ei)^R2{92)\g^2) 



(20) 



with 9i and 62 the angles defining the orientation of the 
observable Ri. We then construct the polynomial TO2 by 
identifying each term appearing in the first of Eqs. (jl4p 
with the correlation function obtained by using two dif- 
ferent sets of angles, [Ox, 62) and {0'i,9'2)- We can then 
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construct the CHSH function as 



S2 



Classically this function is bounded by 1, while there are 
quantum mechanical states that violate this bound by 
up to a factor \/2. In Fig [3] (a) we see how a typical 
CHSH function |iS2| behaves against the amplitude of the 
transverse field B, for a set value of J. For nearly van- 
ishing values of B we have a maximally entangled state 
which maximally violates the bound set by local realis- 
tic theories. As the amplitude of B increases, we see ^2 
asymptotically approaching the local realistic bound of 
1. 

Despite the apparent simplicity of the case studied 
above, the results achieved in the two-spin ring serve as 
useful milestones in the analysis to follow. In fact, we 
can extend our investigation to the genuinely multipar- 
tite scenario by using the information on the ground state 
of the Ising ring under scrutiny that can be gathered by 
means of the approach described in Sec. iHland applying 
it to the formalism for multipartite Bell inequalities given 
in Sec. IIIII Our approach goes as follows: although cer- 
tainly possible, the calculation of the correlation function 
in terms of the fermionic operators putting H into the 
form Hg given in Eq. ([6|) is not straightforward, especially 
under the scalability viewpoint. Indeed, as significantly 
remarked by Osborne and Nielsen [l^, calculating cor- 
relation functions for more than two particles is, in gen- 
eral, a daunting problem. One should, in fact, work out 
the expression of Ri{0i) in the fermionic representation, 
keep in consideration the anti-commutation rules obeyed 
by fermions and re-cast the tensor product of A'^-particle 
local observables in terms of the operators bk and b\. Ex- 
pectation values over the ground state |gs^) of an TV-spin 
ring could be finally evaluated. Instead, we decided to 
use the information on the energy A jv obtained by means 
of the apparatus set up in Sec. HIlso as to determine the 
form of the ground state of the system as expressed in 
terms of the elements of the computational basis. This 
is easily and computationally very conveniently done by 
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FIG. 2: (Color online). Comparison of the A''-concurrence NC 
corresponding to the ground state of a ring of TV = 2, 4 and 
6 spins (curves from top to bottom one) at zero temperature. 
The A'-concurrence is plotted against B and J. 




FIG. 3: (Color onUne). (a) Multipartite non-locality of the 
ground state of a ring of an even number of spins against the 
amplitude of the magnetic field B and at J = 5. The solid 
line is for A'' = 6, the dashed line for A'^ = 4 and, finally, the 
dotted one is for N — 2. (b) Same as in panel (a) but for an 
odd number of spins. The dotted line is for A = 3, the dashed 
for N — 5 while the squares show the results for A = 7. 



solving the set of equations 

{H -ANt2«) \gsN) = 0, 



(22) 



where Ajv is now a known quantity, 1 is the identity ma- 
trix and is the identically null vector. Once the struc- 
ture of the ground state is determined, the evaluation of 
the necessary correlation functions is largely simplified. 
To the best of our knowledge, our study is the first to 
address multi-spin correlations over the ground state if 
an Ising ring. We have thus studied the behavior of the 
Svetlichny function against the amplitude of the mag- 
netic field B at a set value of J (we remind that the 
important parameter, in this problem, is the ratio J / B). 
In Figs. 13] (a) and (b) we show the results of our investi- 
gation for rings of up to 7 spins. It is worth stressing that 
larger rings can well be studied by our method. Given 
the determination of the ground state has been reduced 
to the solution of the eigen-equation Eq. (1^ we require 
only the solution of 2^^ linear equations. It is in the eval- 
uation of the multipoint correlation functions that care 
must be taken. Here the number of correlation functions 
to be evaluated scales with 2^, and we must then op- 
timize the generalized Svetlichny polynomials Eq. (jlSp 
with respect to 2A parameters. This leads to escalating 
computational cost making rings larger than A^ = 7 in- 
tractable. We have split the behavior of the 5„ functions 
against the parity of the length of a ring, thus defining 
Seven and Sodd- As it will be shown later, parity appears 
to be relevant in the investigation of the features here 
under scrutiny. 

The dependence of the the two parity-related 
Svetlichny functions on B is mutually fairly similar and 
overall consistent with the behavior revealed for 82'- at 
large values of J/B, the ground state of the ring exhibits 
strong non-locality, up to the maximum allowed value 
set by the number of spins being considered. As B grows 
(for a set value of J), the non-local nature of |gs„) de- 
teriorates. Remarkably, though, while the even-parity 
case approaches the local realistic bound only asymptot- 
ically with B, a finite value of its amplitude is sufficient 
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to set \Sodd\ = 1 exactly. Although it is difficult to es- 
tablish an explicit relation at this level, we believe that 
such a parity-induced effect is related to the differences 
in the way the ground state of an Ising ring has to be 
considered in the even- and odd-number of spins cases, 
as highlighted in Sec. IIIII when the ring Hamiltonian was 
diagonalized. However, it should be remarked that longer 
chains of an odd number of spins require a larger B to 
reach the local realistic bound, while in a long chain, the 
addition of a single spin would make a very small differ- 
ence, as it is intuitively reasonable to expect. Therefore, 
one can extrapolate to the thermodynamical limit the 
behavior of non-locality and claim that the parity effect 
highlighted above would eventually disappear. As a side 
remark, the reduced state of any two-spin subsystem, re- 
gardless of the length of a ring, despite being entangled 
do not violate Eq. (l2Tt . 



V. BEHAVIOR OF NON CLASSICAL 
CORRELATIONS: THERMAL EFFECTS 

We now investigate the behavior of non-classical corre- 
lations in the non-zero temperature case. Using Eq. ([3]) 
we can determine the thermal state of the Hamiltonian. 
In general, this necessitates of the full spectrum of eigen- 
states the Hamiltonian, whose determination for a large 
number of spins becomes a very difficult task to achieve. 
However, very useful information comes from the study 
of the energy spectrum of H. Fig. 2] shows the spec- 
trum of a ring of two, three and four spins for J = b 
and against increasing field B [30[. It is seen that the 
splitting between the two lowest energy states increases 
with the strength of the magnetic field. In the region of 
B <^ J , the ground and first excited energy levels are 
almost degenerate and a much larger energy gap sepa- 
rates the quasi-degenerate doublet from the rest of the 
spectrum. This implies that an excellent approximation 
of the thermal state of the ring comes from considering 
an admixture of |gs^) and jfirstTv) —b]. Igs^v) with coeffi- 
cients determined by the Boltzmann factor corresponding 
to the energies Ajv and A]y= (first at] iJff jfirstAr) such as 

/ /3Ajv <9a1, \ 

e(r)~ ( e-— |gs^) (gsjvl -(-e-— IfirstAr) (firstwl j /Z. 

(23) 



(a) (b) (c) 




FIG. 4: (Color online). Energy spectrums for a spin ring of 
N = 2 [panel (a)], 3 [panel (b)] and 4 [panel (c)] against B. 
We have taken J = 5. 



Clearly, the approximation is legitimate as long as 
the state (|23p remains quasi-normalized, that is 

(exp[-^] + exp[-^])/Z ~ 1. The proximity of the 
first two energy levels for small values of B implies that 
even a small amount of thermal energy is enough to fully 
mix the two lowest eigenstates, which might well result 
in the wash-out of the quantum correlations within the 
system. On the other hand, for increasing values of -B, a 
larger temperature is required to mix the eigenstates. 

We have examined the entanglement properties for 
even values of N and N = 3 (which, we remark, is 
the only odd case treatable through tripartite negativity) 
against the temperature. For B <^ J, entanglement sud- 
denly disappears at finite values of T. On the other hand, 
for B K J, quantum correlations appear to cope much 
better against the effects of a non-zero temperature: as 
B gets larger, the state of the ring is able to withstand 
larger values of T before we starting to decay. Yet, this 
effect is accompanied by a reduction in the amount of 
entanglement shared by spins. Again, in order to by- 
pass the problems inherent in entanglement measures for 
mixed multipartite states, we use the tools provided by 
multipartite non-locality tests. Figs. [5] (a)- (c) show the 
behavior of the generalized Svetlichny functions for a few 
significant cases. Evidently, as B grows, while paying in 
terms of degree of violation of the appropriate Bell-like 
inequality, we significantly enlarge the range of T within 
which the correlations within the state of a ring cannot 
be described by a local realistic theory. 

This protection effect can be understood as follows: 
the increasing magnetic field tends to align the spins 
in the ring, thus affecting the strength of the spin-spin 
correlations. On the other hand, the tendency to mu- 
tual alignment of the spins makes them less sensitive 
to the fluctuations induced by thermal effects. At the 
same time, a growing B makes the spectral gap between 
ground and first excited state of the ring wider, thus 
making the admixture between the ground state and the 
rest of the spectrum less likely to occur. We can link 
the effects analyzed above to the magnetization of the 
ring so as to provide a much neater connection between 
the correlation properties of the system and its intrinsic 
macroscopic features, which is certainly a very interest- 
ing point to assess. Let us concentrate on the two-spin 
ring case, which is non-trivial enough to encompass the 
general features we are interested in. We consider the 
expectation value of the global magnetization operator 
M.2= Y^j=i O'^s^ the state of a two-spin ring at tem- 
perature T. Clearly, the modulus of this function grows 
monotonically from zero to one at a rate slowed down 
by the increasing temperature. Here, we are interested 
in the connection between non-classical correlations and 
the global magnetization {M.2)- Due to the dimension 
of the problem at hand, we can equally use negativity 
or the Bell-like function 52 and we focus on the former 
simply for easiness of calculations. In Fig. [5] we plot N2 
against {M.2) for a range of values of the temperature. 
At T = 0, small values of the magnetization (correspond- 



7 



(a) (b) (c) 




FIG. 5: (Color online). Behavior of Svetlichny functions Sn for the state of a ring with J = 5 and T ^ 0. We address the case 
of = 2, 3 and 4 (panel (a), (b) and (c) respectively) as representative of the general trend exhibited by arbitrary A'' values. 
Each curve in the same plot is for for a specific choice of B £ [1, 10]. Regardless of N, as B increases the degree of violation of 
the appropriate Bell-like inequality decreases. However, the resilience to spoiling effects due to temperature increases. 



ing to large values of the ratio J/B) correspond to a 
very large degree of entanglement, in agreement with the 
analysis so far [ll| . A growing magnetization sees A2 de- 
creasing to zero. The behavior changes abruptly for even 
small non-zero values of T, when we start from an ini- 
tially maximally entangled state and as we increase the 
magnetization of the system we see a decrease in non- 
classical correlations. For the bipartite case we see the 
correlations become purely classical only when the sys- 
tem is fully saturated by the magnetization, i.e. when 
M = 1 or equivalently B 00. Qualitatively the results 
hold in the higher dimensional cases. While the study of 
this figure against temperature reinforces the features re- 
vealed above (entanglement decays more slowly at larger 
values of B, although the maximum achievable degree of 
negativity decreases) , we can also infer that larger values 
of the magnetic field are required in order to optimize the 
entanglement within the ring, if T grows. This is a clear 
manifestation of the correlation-by- alignment effect dis- 
cussed before: the thermal randomness is counteracted 




FIG. 6: (Color online). Two-spin negativity against magneti- 
zation and temperature for J = 5. We have taken T = 1, .., 20 
(increasing at unit steps). The entanglement for magnetiza- 
tion going toward 1 stays at finite non-zero values only be- 
cause we have considered B G [0, 100]. By enlarging this 
range, smaller final values of negativity would be reached. In 
fact, we have A/i ^ as {M2) 1. 



by the order induced by a strong magnetic field. How- 
ever, if the alignment induced by B exceeds the threshold 
determined by the degree of thermal disorder, quantum 
correlations start to vanish as the spins tend to be only 
classically correlated (as expected in an Ising ring beyond 
the critical point). A qualitatively similar analysis holds 
for longer rings. By controlling the macroscopic features 
of magnetization, one could thus determine the degree of 
quantum correlations set among the elements of the ring. 

VI. CONCLUSIONS 

We have studied the non-classical correlations natu- 
rally occurring in the ground and thermal state of an 
Ising ring in transverse magnetic field. The focus of our 
analysis was the behavior of multipartite entanglement 
and non-locality as quantified by the violation of ad hoc 
many-spin Bell- like inequalities [28| . At zero tempera- 
ture, while even-length rings cannot be described by lo- 
cal realistic theories even at the limit of asymptotically 
infinite strength of the magnetic field, odd-length ones 
behave differently: there is always a finite choice of B 
beyond which the bound imposed by local realism is sat- 
isfied. We have been able to exactly tackle multipar- 
tite correlation functions in such a many-body problem, 
therefore going significantly far from the standard anal- 
ysis of the behavior of quantum correlations in interact- 
ing many-body models. At non-zero temperature, we 
have found that the thermal randomness that would de- 
stroy the non-classical correlations found in the ground 
state of the model can be counteracted simply by set- 
ting a large enough magnetization of the ring: a large 
magnetic field, in fact, while decreasing the degree of vi- 
olation of Bell-like inequalities, will make the generalized 
Svetlichny functions used for non-locality test resilient 
against thermal spoiling effects. This suggests that a gen- 
uinely multipartite quantum-correlated state in an Ising 
ring can be generated simply by controlling the macro- 
scopic properties of the system, i.e. magnetization and 
temperature. We believe our results will trigger further 
interest in examining entanglement and non-locality in 



8 



many-body interacting systems. able discussions. We acknowledge financial support from 

DEL and the UK EPSRC (EP/G004579/1), 

Acknowledgments 

We thank Dr. Nicolas Brunner, Dr. Tony ApoUaro, 
Prof. Sougato Bose and Prof. Vlatko Vedral for invalu- 



[1] S. Bose Phys. Rev. Lett. 91, 207901 (2003). 

[2] M. C. Arnesen, S. Bose, and V. Vedral, Phys. Rev. Lett. 

87, 017901 (2001). 
[3] A. Bayat, P. Sodano, and S. Bose, Phys. Rev. B 81, 

064429 (2010). 
[4] S. Bose, Contemporary Physics 48, 13 (2007). 
[5] R. Horodecki, P. Horodecki, M. Horodecki and K. 

Horodecki, Rev. Mod. Phys. 81, 865 (2009). 
[6] S. Popescu, Nature Phys.' 6, 151 - 153 (2010). 
[7] M. Ansmann et al, Nature 461, 504-506 (2009). 
[8] N. Brunner and P. Skrzypczyk, Phys. Rev. Lett. 102, 

160403 (2009); J. Allcock, N. Brunner, N. Linden, S. 

Popescu, P. Skrzypczyk and T. Vertesi Phys. Rev. A 80, 

062107 (2009). 
[9] L. Viola and S. Lloyd, Phys. Rev. A 58, 2733 (1998). 
[10] J. G. Oliveira, R. Rossi, and M. C. Nemes Phys. Rev. A 

78 044301 (2008). 
[11] P. Stelmachovic and V. Buzek, Phys. Rev. A 70, 032313 

(2004). 

[12] D. Gunlycke, S. Bose, V. Kendon and V. Vedral, Phys. 
Rev. A 64, 042302 (2001). 

[13] A. Montakhab and A. Asadian, e-print arXiv: 10 09. 13571 

[14] D. M. Greenberger, M.A. Home and A. Zeilinger, Bell 
Theorem, Quantum Theory, and Conceptions of the Uni- 
verse Ed. M. Kafatos (1989 Dordrecht: Kluwer) pg. 69. 

[15] E. Leib, T. Schuhz, and D. Mattis, Ann. Physics 16, 407 
(1961). 

[16] T. J. Osborne and M. A. Nielsen, Phys. Rev. A 66, 

032110 (2002). 
[17] P. Jordan and E. Wigner, Z. Physik 47 631 (1928). 



[18] E. Rico Ortega, Quantum correlations in (1 + 1)- 

dimensional systems, Ph.D. Thesis, Universidad de 

Barcelona (2008). 
[19] Peres, Phys. Rev. Lett. 77, 1413 (1996); M. Horodecki, 

P. Horodecki, and R. Horodecki, Phys. Lett. A 223, 1 

(1996). 

[20] C. H. Bennett, D. P. DiVincenzo, J. A. Smolin, and W. K. 

Wooters, Phys. Rev. A 54, 3824 (1996); W. K. Wootters, 

Phys. Rev. Lett. 80, 2245 (1998). 
[21] C. Sabm and G. Garci'a-Alcafne, Eur. Phys. J. D 48, 

435-442 (2008). 
[22] A. Wong and N. Christensen, Phys. Rev. A 63, 044301 

(2001). 

[23] G. K. Brennen and S. S. Bullock, Phys. Rev. A 70, 
052303 (2004). 

[24] J. S. BeU Physics 1, 195 (1964); J. F. Clauser, M. A. 
Home, A. Shimony, and R. A. Holt, Phys. Rev. Lett. 23 
880 (1969). 

[25] J. Cereceda, Phys. Rev. A 66, 024102 (2002). 
[26] G. Svetlichny Phys. Rev. D. 35, 3066 (1987). 
[27] N. D. Mermin, Phys. Rev. Lett. 65, 1838 (1990); D. N. 

Klyshko, Phys. Lett. A 172, 399 (1993). 
[28] D. GoUins, N. Gisin, S. Popescu, D. Roberts, and V. 

Scarani, Phys. Rev. Lett. 88 170405 (2002). 
[29] M. Seevinck and G. Svetlichny, Phys. Rev. Lett. 89, 

060401 (2002). 

[30] Analogous spectra hold for larger rings, although the 
quickly growing number of energy levels make their visu- 
alization difScult. 



